Next Article in Journal
Exploring the Origin of Maximum Entropy States Relevant to Resonant Modes in Modern Chladni Plates
Next Article in Special Issue
Effects of Vaccination Efficacy on Wealth Distribution in Kinetic Epidemic Models
Previous Article in Journal
Power Graphs of Finite Groups Determined by Hosoya Properties
Previous Article in Special Issue
On the Mortality of Companies
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Analysis of Individual High-Frequency Traders’ Buy–Sell Order Strategy Based on Multivariate Hawkes Process

1
Department of Mathematical and Computing Science, School of Computing, Tokyo Institute of Technology, Yokohama 226-8502, Japan
2
Institute of Innovative Research, Tokyo Institute of Technology, Yokohama 226-8502, Japan
3
Sony Computer Science Laboratories, Tokyo 141-0022, Japan
*
Author to whom correspondence should be addressed.
Entropy 2022, 24(2), 214; https://doi.org/10.3390/e24020214
Submission received: 7 January 2022 / Revised: 26 January 2022 / Accepted: 28 January 2022 / Published: 29 January 2022
(This article belongs to the Special Issue Three Risky Decades: A Time for Econophysics?)

Abstract

:
Traders who instantly react to changes in the financial market and place orders in milliseconds are called high-frequency traders (HFTs). HFTs have recently become more prevalent and attracting attention in the study of market microstructures. In this study, we used data to track the order history of individual HFTs in the USD/JPY forex market to reveal how individual HFTs interact with the order book and what strategies they use to place their limit orders. Specifically, we introduced an 8-dimensional multivariate Hawkes process that included the excitations due to the occurrence of limit orders, cancel orders, and executions in the order book change, and performed maximum likelihood estimations of the limit order processes for 134 HFTs. As a result, we found that the limit order generation processes of 104 of the 134 HFTs were modeled by a multivariate Hawkes process. In this analysis of the EBS market, the HFTs whose strategies were modeled by the Hawkes process were categorized into three groups according to their excitation mechanisms: (1) those excited by executions; (2) those that were excited by the occurrences or cancellations of limit orders; and (3) those that were excited by their own orders.

1. Introduction

To gain a deeper understanding of the mechanisms of financial markets, it is necessary to clarify the order strategies of individual market participants. In financial markets, recent developments in information technology have made it possible to track the transactions of individual market participants in detail. These technological advances have led to the analysis of the trading strategies of individual market participants and how these strategies affect financial markets. For example, Odean [1], and Grinblatt and Keloharju [2], reported the relationship between historical returns and market participants’ decisions to buy and sell stocks. The position management strategies of individual market participants were analyzed based on the data, which confirmed that these strategies actually affected market prices in the near future [3]. Individual strategies for placing buy and sell orders in response to market price changes were approximated using a simple mathematical model, and the basic statistical properties of financial Brownian motion were theoretically derived based on the kinetic theory in a manner parallel to traditional statistical physics [4,5,6].
In particular, high-frequency traders (HFTs) have recently attracted attention. HFTs are algorithmic traders who can react to market changes in milliseconds and place or cancel, buy and sell orders at high frequencies [7]. Because of the development of information technology, they have a large presence in financial markets around the world. In fact, HFTs accounted for 68.3% of the total trading volume in the stock market [8]. Furthermore, HFTs currently account for the majority of orders shown in the order book [9,10,11]. The availability of high-frequency trading data has triggered the academic study of HFTs [12,13]. Previous studies have generally agreed that HFTs make market spreads smaller and enhance market liquidity [8,14,15]. As an indicator for predicting the short-term volatility of market prices from the order book information, the volume-synchronized probability of informed trading (VPIN) has been proposed and actively studied [16,17,18,19,20]. In addition, informed trading using the advantage of information such as public news and confidential information has been studied using high-frequency data [21,22,23,24,25]. We believe that it is crucial to gain a deep understanding of their trading behavior in current financial markets, where HFTs provide most of the liquidity.
In this study, we used a multivariate Hawkes process to investigate the processes used by individual HFTs for generating sell and buy limit orders in the USD/JPY forex market, and clarified when each HFT placed buy–sell limit orders. The Hawkes process is a type of non-homogeneous Poisson process proposed by Hawkes [26]. As will be explained later, it is characterized by an intensity function which determines the probability of the occurrence of an event in a point process. It utilizes an excitation term that is affected by past events, and can describe a point process associated with past events. Similar ideas have been independently introduced for financial markets to explain the strong correlation to past events, such as the “autoregressive conditional duration model” [27] and “self-modulation processes” [28]. The Hawkes process is a useful model for interpreting financial phenomena, in which many factors interact to produce complex aspects. In this paper, we show that it is also useful for interpreting the behavior of HFTs. Specifically, we introduced a multivariate Hawkes process in which the process of generating HFTs’ buy–sell limit orders is mutually excited by a total of eight events, such as the creation of limit orders, the cancellation of limit orders, and execution in the order book, showing that the order behaviors of many HFTs can be modeled by the Hawkes process.
Hawkes processes [29] have various applications in the financial field, such as those related to volatility clustering [30], market activity and risk [31,32,33], and market impact [34]. In particular, the Hawkes process has been actively employed as an approach to the dynamic description of order books, where a set of order types is specified and a multivariate Hawkes process is fitted to their timestamps [35,36,37,38,39,40,41]. However, there has been no study that used a multivariate Hawkes process to investigate the order generation processes of individual HFTs. In today’s financial markets, where the majority of order books are made up of HFTs’ orders, our empirical results provide new information from a more microscopic perspective. We believe that this study shed light on how HFTs provide liquidity to the market.
The remainder of this paper is organized as follows. Section 2 explains the datasets and describes the HFTs that were analyzed in this research. Section 3 introduces the multivariate Hawkes process and describes the method used for parameter estimation. In Section 4, a clustering analysis of 134 HFTs is introduced to categorize their strategies based on the estimated Hawkes’ parameters. In Section 5, we discuss our results.

2. Data

First, we provide a basic description of the order data for the USD/JPY forex market (EBS market), along with individual trader IDs (see Section 2.1). We then define the HFTs in this market (see Section 2.2) and show some examples to explain how their buy–sell limit order generation is linked to changes in the order book (see Section 2.3).

2.1. EBS Market Data Description

In this study, we used high-frequency data for the USD/JPY forex market provided by the EBS. EBS is an interbank forex market and one of the largest financial platforms in the world. Because it is an interbank market, most market participants are professional traders from banks, hedge funds, and other financial institutions, and our forex dataset contains their trading data. Our dataset contains information from five days (from 21:00 GMT on 5 June 2016 to 21:00 GMT on 10 June 2016), with a total of approximately 2.8 million orders and a transaction volume of USD 68 billion corresponding to this period. Table 1 shows an example of the raw data we used. The data for each of the 2.8 million orders contained not only the order type, price, volume, and timestamp (in milliseconds), but also an anonymized trader ID that could identify who submitted the order. Using these trader IDs, we could track individual traders’ full orders in milliseconds. In addition, the minimum price unit that a trader could submit was JPY 0.005, and the minimum transmission volume was USD 1 million.
The EBS market is open 24 h a day from Monday morning to Friday at midnight, and trading is conducted via a double auction system in the order book. Figure 1 shows a schematic of the trading in the order book, where the horizontal axis is the price and the vertical axis is the volume. There are six order types for trading: buy/sell limit orders, buy/sell cancel orders, and buy/sell market orders. Limit orders are submitted at the trader’s desired price and remain in the order book until traded or cancelled. Cancel orders are submitted by a trader to cancel a limit order that they previously submitted. Market orders are submitted at the current best limit price. Transactions that are executed by buy market orders are called hit sell transactions, and transactions executed by sell market orders are called hit buy transactions (see Figure 1). If the best price worsens (e.g., the best sell limit price becomes higher) before the market receives the market order at the best price, the market order is automatically invalidated. In fact, this study found that 79.5% of the market orders were invalidated without being executed.
Figure 2a shows the average trading price per 10 min window over the 5 days we analyzed. During this period, there are no market crashes or spikes. Figure 2b shows the number of each type of order per day, which looks stable.

2.2. Definition of HFTs

“HFTs” is a general term for traders who place and cancel orders at high speed and high frequency according to an algorithm, but there are various definitions. In this study, we define an HFT as a trader who places both buy and sell limit orders and presents an average of 500 or more limit orders per day following the previous report written by a researcher from EBS [9]. Based on this definition, the number of HFTs was 134 out of the 1031 traders included in this 5-day data set. These 134 HFTs accounted for 89.6% of the market’s total number of limit orders.
Figure 3a shows the histogram of the minimum time interval between orders for each HFT. There is no description in the data to identify whether the ID is a human or a computer; however, Figure 3a shows that most of the intervals are within 0.1 s, which are difficult for a human to execute.
In Figure 3b, we plot the number of HFTs and non-HFTs participating in the market every hour, indicating that the number of HFTs is relatively stable compared to non-HFTs. Figure 3c shows the percentage of limit orders placed by HFTs every hour, demonstrating that the majority of the limit orders are provided by the HFTs.

2.3. Basic Properties of HFTs

In this study, we focused on the limit order generation process of HFTs, which accounted for the majority of limit orders in the order book. Naturally, the order strategies (i.e., the processes used to submit limit orders) of HFTs differed from every algorithm. However, it is natural for them to see the quotes in the order book when submitting their limit orders. Figure 4a,b plot the numbers of buy–sell limit orders per 10 min window for three HFTs, respectively, and Figure 4c plots the numbers for six types of orders in the order book. From Figure 4, we can observe that the numbers of buy–sell limit orders from the three HFTs increased or decreased simultaneously and tended to be in sync. More interestingly, the numbers of these HFTs’ buy–sell limit orders tend to be in sync with the numbers for each type of order in the order book where all market participants’ orders are submitted. Since the above synchronization phenomenon was confirmed for many HFTs, we believe that many HFTs react instantaneously to some changes in the market when submitting limit orders.

3. Method

The preceding section showed that the limit order generation processes of the HFTs tended to be in sync with traders’ orders or orders in the order book. To clarify how 134 HFTs’ buy–sell limit orders react to the other order events, their order processes are modeled by the multivariate Hawkes process. In this section, we introduce the multivariate Hawkes process that we used in this study (see Section 3.1) and explain the parameter estimation method based on maximum likelihood estimation (see Section 3.2). We then describe the validity of the estimation results (see Section 3.3).

3.1. Model

This section presents an overview of the Hawkes process and introduces our model.

3.1.1. Mathematical Notation

Let us consider point process { t i } , which is a sequence of non-negative random variables such that i N , t i < t i + 1 . For point process { t i } , the conditional intensity function is defined as follows [42]:
λ t H t = lim Δ t 0 P N ( t + Δ t ) N ( t ) = 1 H t Δ t
where P { A | B } represents the probability of A under condition B, N ( t ) is the cumulative number of event occurrences at time t (i.e., a counting process), and H t is the history of the process up to time t containing a sequence of event times { t i } (i.e., a filtration). As can be seen from the definition, λ t H t Δ t represents the probability of an event occurring in time interval [ t , t + Δ t ) . Here, we use the shorthand notation λ ( t ) λ ( t H t ) , assuming the history up to time t, and we call the Poisson parameter λ ( t ) an intensity function.

3.1.2. Overview of Hawkes Process

The Hawkes process is a point process in which the intensity function is affected by the occurrence of past events. Let { t i } be a point process, and N ( t ) be the associated counting process, and the intensity function of the generalized Hawkes process is defined as follows [43]:
λ ( t ) = c + t ϕ ( t s ) d N ( s ) = c + t i < t ϕ t t i
where c is a positive constant showing a base intensity, and ϕ ( t ) is a kernel function that expresses the effect of event t i from the past on the current intensity [44]. There are various types of kernel functions, and their properties have been well studied. In this study, we applied an exponential kernel α e β t , which is a popular kernel function that was originally proposed by Hawkes [26]. We call this Hawkes process a univariate Hawkes process because it is affected by its own events.
The Hawkes process can be extended to a multivariate model in which several types of point processes interact with each other. Let { t i } { { t 1 , i } , { t 2 , i } , , { t M , i } } be M-variable point processes, and N ( t ) = { N 1 ( t ) , N 2 ( t ) , , N M ( t ) } be the associated counting process, the intensity function of a multivariate Hawkes process for point process { t n , i } is defined as follows [43]:
λ n ( t ) = c n + m = 1 M t ϕ n , m ( t s ) d N m ( s ) = c n + m = 1 M t m , i < t ϕ n , m t t m , i
As in the case of the univariate Hawkes process, c n is a positive constant showing a base intensity, and ϕ n , m ( t ) is a kernel function that expresses the effect of event t m , i from the past on the current intensity. In this case, ϕ n , m ( t ) = α n , m e β n , m t , with positive constant parameters α n , m and β n , m , and the intensity function for point process { t n , i } is given as follows:
λ n ( t ) = c n + m = 1 M t α n , m e β n , m ( t s ) d N m ( s ) = c n + m = 1 M t m , i < t α n , m exp ( β n , m t t m , i )
Figure 5 is a schematic of this intensity function (Equation (4)). The intensity function, λ n ( t ) , increases α n m at event time t n , i and exponentially decays with a time constant of 1 / β n m . Thus, λ n ( t ) is excited not only by its own events, but also by other events, and the multivariate Hawkes process can represent such mutual interactions.
The quantity, ρ n , m , expressed in the following equation (Equation (5)) in the case of an exponential kernel is called the branching ratio [45,46]. This is the expectation of the number of occurrences of event n caused by the occurrence of event m. A larger value for this number represents a greater impact of event m on event n, and this value is an important quantity for interpreting the Hawkes process:
ρ n , m α n , m β n , m = t m , i α n , m exp β n , m t t m , i d t

3.1.3. Trader Model

In this study, the buy and sell limit order processes of 134 HFTs are modeled by eight-variable Hawkes processes with exponential kernels that are excited by a total of eight-point processes. These eight types were the target HFTs’ sell limit (TS) and buy limit (TB), and six types of orders in the order book: sell limit (SL); buy limit (BL); sell cancel (SC); buy cancel (BC); hit sell (HS); and hit buy (HB).
Let { t i } { { t T S , i } , { t T B , i } , { t S L , i } , { t B L , i } , { t S C , i } , { t B C , i } , { t H S , i } , { t H B , i } } be an eight-variable point process, the intensity functions for { t T S , i } and { t T B , i } are given as follows:
λ T S ( t ) = c T S + m M t m , i < t α T S , m exp ( β T S , m t t m , i )
λ T B ( t ) = c T B + m M t m , i < t α T B , m exp ( β T B , m t t m , i )
where M { T S , T B , S L , B L , S C , B C , H S , H B } . We assume that the above intensity functions (Equations (6) and (7)) represent the buy and sell limit order processes of each HFT and examine which events affect their order generation processes. Hereafter, the abbreviations listed in Table 2 are used for the order events.

3.2. Parameter Estimation Using Maximum Likelihood Estimation

For Equations (6) and (7), we apply the maximum likelihood estimation method for the parameter estimation. Because the functional forms of the intensity functions are the same for { t T S , i } and { t T B , i } , we solve the log-likelihood functions in the following manner. For the point process { t T S , i } , the log-likelihood function in time interval [ 0 , T ] is given by the following [47]:
log L ( c T S , α T S , β T S ) = 0 T λ T S ( t ) d t + i = 1 n log λ T S t i = c T S T + m M α T S , m β T S , m t m , i < T { e x p ( β T S , m ( T t m , i ) 1 } + i = 1 n log c T S + m M t m , j < t T S , i α T S , m e x p ( β T S , m ( t T S , i t m , j ) )
The same formulation is also applied for log L ( c T B , α T B , β T B ) .
These log-likelihood functions are differentiable by each parameter, and we optimized them by Adam [48], which is a type of gradient descent method, to obtain the maximum likelihood estimators. Here, the initial values are ( c T S , α T S , β T S ) = ( c T B , α T B , β T B ) = ( 0.1 , 0 . 1 , 10 ) , and the various parameters required for Adam are set following the values in the original reference [48]. Here, α T S , β T S , α T B , and β T B are vectors of eight variables (e.g., α T S = ( α T S , T S , α T S , T B , , α T S , H S ) ).
It is also known that the computation of the gradients of the log-likelihood function of the Hawkes process usually requires the computation of O ( N 2 ) . However, using a recursive formulation that can be used when the kernel function is an exponential function, we perform maximum likelihood estimation with O ( N ) computational complexity (see Ogata [49] for the recursive formulation).
In addition, the HFTs do not always continuously place orders. Therefore, if an HFT did not place any orders for more than 15 min, we considered such period as not participating in the trade and performed the maximum likelihood estimation for each trader by ignoring such inactive periods.

3.3. Validity of Estimation Results

Based on a residual analysis [50], we assessed the goodness-of-fit of the point process model. Let the intensity function for point process { t i } be λ ( t ) , then the sequence, { τ i } , of random variables, where each element is transformed by τ i 0 t i λ ( t ) d t , has the distribution of a stationary Poisson process with intensity 1, and the transformed residual Δ τ i τ i + 1 τ i has an exponential distribution with the unit mean.
Therefore, if the estimated HFT intensities, λ ^ T S ( t ) and λ ^ T B ( t ) , are good approximations of the true intensities, λ T S ( t ) and λ T B ( t ) , respectively, then the transformed residuals, Δ τ ^ T S , i and Δ τ ^ T B , i , are expected to follow the exponential distributions with the unit mean. Here, the transformed residuals are defined as Δ τ ^ T S , i t T S , i t T S , i + 1 λ ^ T S ( t ) d t and Δ τ ^ T B , i t T B , i t T B , i + 1 λ ^ T B ( t ) d t , respectively, which can be derived in the case of Δ τ ^ T S , i as an example (Equation (9)) because the intensity function has the same form as described above:
Δ τ ^ T S , i = c ^ T S Δ t T S , i m M t m , i < t T S , i + 1 α ^ n , m β ^ n , m ( e x p ( β ^ T S , m ( t T S , i + 1 t m , i ) ) e x p ( β ^ T S , m ( t T S , i t m , i ) ) )
Figure 6a shows the cumulative distribution of the three HFTs’ original residuals, Δ t T S , i , and Figure 6b shows the cumulative distribution of the same three HFTs’ transformed residuals, Δ τ ^ T S , i . We confirmed that the transformed residuals for the three HFTs approximately followed the exponential distribution with the unit mean, which implied that the intensities of { t T S , i } were properly approximated. Because not all of the HFTs’ order generation processes could be modeled by the Hawkes process proposed here, the above operations were performed for the sell and buy order processes of the 134 HFTs.
Here, we apply the Kullback–Leibler divergence between the distribution of the transformed residuals and the exponential distribution with the unit mean as a criterion to determine whether the approximations of the intensities of each HFTs’ { t T S , i } and { t T B , i } are appropriate. Since we estimate the intensity functions of the two-point processes for each trader, we calculate the sum of the respective Kullback–Leibler divergences ( D T S + T B K L ), as defined by Equation (10) below:
D T S + T B K L = j p j Δ τ ^ T S , i log p j Δ τ ^ T S , i q j Δ τ ^ T S , i + j p j Δ τ ^ T B , i log p j Δ τ ^ T B , i q j Δ τ ^ T B , i
Here, q j ( τ ) is the discrete exponential distribution with the unit mean, and p j τ is the sampling distribution. The bin size of the discrete distribution is assumed to be 1. The threshold of accepting D T S + T B K L error will be determined in the next section.

4. Results

In this section, we first show that the Hawkes process introduced here successfully approximated the intensity of the limit order generation process for 104 of the 134 HFTs by applying the method described in Section 3.3 (see Section 4.1). We then categorize the order generation processes of the 104 successfully estimated HFTs into three groups according to their excitation mechanisms, and explain how each group of HFTs places their orders (see Section 4.2).

4.1. D T S + T B K L Calculation Results for All HFTs

Figure 7 shows a histogram of the D T S + T B K L values for the 134 HFTs. From the histogram, there are many HFTs whose values of D T S + T B K L are very close to 0 and the plots are scattered for D T S + T B K L > 0.05 , so we set the threshold as 0.05. The 104 HFTs who fell into the range of D T S + T B K L < 0.05 were considered to be traders whose order generation processes were well modeled, while the rest of the 30 HFTs who fell into the range of D T S + T B K L 0.05 were considered to be traders whose order generation processes were poorly modeled by the Hawkes process introduced here.
With the estimated parameters, the Hawkes process could be simulated using the thinning method [51]. For example, Figure 8 compares the time series of the number of orders per 10 min window for the real data of an HFT with D T S + T B K L = 0.0238 and a simulated time series. It can be confirmed that both sell limit orders (upper figure) and buy limit orders (lower figure) successfully reproduce the behavior of the real data.
On the other hand, Figure 9 shows a comparison of the simulated and real data for an HFT with D T S + T B K L = 0.211 , which was judged not to be properly estimated by the Hawkes process. It can be seen that the deviations from the real data for both sell limit orders (upper figure) and buy limit orders (lower figure) are larger than in the case of Figure 8. The Hawkes process introduced here did not adequately explain the order generation process for this trader with a large error, D T S + T B K L . Because some traders could not be modeled by the Hawkes process, in the following, we report the results of our clustering analysis of the order generation processes of 104 HFTs after excluding 30 traders.

4.2. Results of Clustering Analysis

Here, we categorize the 104 HFTs whose limit order generation process was properly estimated according to the similarity of their excitation mechanisms (the branching ratio), and describe how each group of HFTs placed buy–sell limit orders and provided liquidity to the market. As defined in Equation (5) in Section 3.3, the branching ratio, ρ n , m , is an absolute value that represents the expectation for the number of occurrences of event n caused by the occurrence of event m. To evaluate the excitations of { t T S , i } and { t T B , i } relative to each HFT, we introduced normalized branching ratios ρ ¯ T S , m and ρ ¯ T B , m , which are defined by the following equations, so that the sum is equal to 1:
ρ ¯ T S , m ρ T S , m i M ρ T S , i
ρ ¯ T B , m ρ T B , m i M ρ T B , i
Because both { t T S , i } and { t T B , i } are 8-variable Hawkes processes, 16 normalized branching rates were defined for each HFT. Figure 10 shows a dendrogram of the hierarchical clustering of the 104 HFTs using these 16 variables. For this hierarchical clustering, we used the Ward method [52] to join clusters in the order of the decreasing sum of squares after joining. The vertical axis in Figure 10 represents the distance between clusters with an increase in the sum of squares when clusters A and B are joined, and is defined by the following equation:
Δ ( A , B ) = i A B x i m A B 2 i A x i m A 2 i B x i m B 2 ( 12 )
where | | d | | denotes the Euclidean distance, and m j is the center of cluster j.
Based on the distance between the clusters, we found it reasonable to categorize the HFTs into three groups with the threshold distance around 3, as shown in Figure 10, and designated them as Group A, Group B, and Group C. There were 77 HFTs in Group A, 12 in Group B, and 15 in Group C. The number of clusters becomes larger for a lower threshold distance, however, we confirmed that properties of any smaller groups are quite similar to one of these three groups in the graphical representation of an interaction network to be discussed in the following.
The remainder of this section explains the order events that excited the HFTs in each group to place buy–sell limit orders based on the estimated Hawkes parameters.

4.2.1. Group A

Group A is comprised of 77 HFTs. The total number of limit orders in the 5 days was approximately 850,000, which accounted for 62.5% of the total number of limit orders in the market. Figure 11 shows the quartiles and means of the 16 normalized branching ratios for these 77 HFTs, where (a) represents ρ ¯ T S , m and (b) represents ρ ¯ T B , m . From Figure 11a, it can be seen that the generation of sell limit orders by the HFTs in Group A was most excited by hit sell, which greatly exceeded the excitation from other events. In contrast, Figure 11b shows that their buy limit order generation was most excited by hit buy, which also greatly exceeded the excitation from other events.
Figure 12 illustrates the network graph of the buy and sell limit orders of the HFTs in Group A, along with all types of orders, using these normalized branching ratios. The size of the directed edges of the network is proportional to the mean value of the normalized branching ratio, e.g., edges directed from HS to TS and from HB to TB represent strong excitations.
In addition, because the kernel function of the Hawkes process is an exponential function, the time constant, which is a measure of the response speed of one excitation at a time, is given by β n , m 1 . The mean values of the estimated time constant for the HFTs in Group A are summarized in Table 3. It is suggested that their reaction speed to an event is approximately 0.1 s, which is reasonable for HFTs who trade at very high speeds.

4.2.2. Group B

Group B is comprised of 12 HFTs. The total number of limit orders in the 5 days was approximately 174,000, which accounted for 12.7% of the total number of limit orders in the market. Figure 13 shows the quartiles and means of the 16 normalized branching ratios for these 12 HFTs. From Figure 13a, we can see that their sell limit order generation was excited by sell limit and cancel buy, and from Figure 13b, we can see that their buy limit order generation was excited by buy limit and cancel sell. Unlike Group A, they did not react to execution events but were excited by the generation and cancellation of limit orders in the order book.
Figure 14 shows the interaction network of the buy and sell limit orders of the HFTs in Group B, along with all types of orders, using the mean of the normalized branching ratios, as in Figure 12.
The mean values of the time constants for each event are summarized in Table 4. As in the case of HFTs in Group A, these values suggest that the reaction speed to events were to be measured in milliseconds.

4.2.3. Group C

Group C is comprised of 15 HFTs. Their total number of limit orders in the 5 days was approximately 95,000, which accounted for 6.9% of the total number of limit orders in the market. From Figure 15a, it can be seen that the sell limit order generation of the HFTs in Group C was most strongly excited by their own buy limit, and was also excited by the sell limit and cancel buy. On the other hand, Figure 15b shows that their buy limit order generation was most strongly excited by their own sell limit, but was also excited by buy limit and cancel buy.
Figure 16 shows the interaction network of the buy–sell limit orders of the HFTs in Group C, along with all types of orders, as in Figure 12 and Figure 14. The HFTs’ sell/buy limit orders interacted with each other.
The mean values of the time constants for each event are summarized in Table 5. The time constants of the excitations from TS to TB and from TB to TS were larger than 10 s, suggesting that the excitations were sustained for a very long time compared to those previously observed.

5. Conclusions

In summary, we introduced a multivariate Hawkes process to model the limit order generation processes of individual HFTs participating in the USD/JPY foreign exchange market for 5 days and analyzed their limit order generation mechanisms. First, we confirmed that an eight-variable Hawkes process, which consisted of each HFTs’ own buy–sell limit orders and the six types of orders in the order book, could adequately model the limit order generation processes of 104 of the 134 HFTs. Then, we categorized the 104 properly modeled HFTs into three categories based on the similarity of the excitation mechanisms measured by the parameter values of the Hawkes process. As a result, we confirmed that the majority of the HFTs in our dataset reacted to the execution of trades, while 12 of the 134 HFTs only reacted to limit orders and 15 of the 134 HFTs reacted to their orders. By evaluating the time constants of the estimated excitations of individual HFTs, we found that many HFTs responded to the most recent change in the order book in a very short time, by placing or canceling new orders. Since HFTs currently account for the majority of limit orders shown in the order book, the results of this analysis provide more microscopic insight into the dynamics of the order book than previous studies.
The following issues will be studied in the future as a generalization of the present work. The first goal is to clarify the limit order generation processes of the remaining 30 HFTs who could not be adequately modeled by the present analysis. The Hawkes process adopted in this study only included the impact of the occurrence of a recent order event and ignored important financial market influences such as the volume of orders, market price fluctuations and trends, and the positions of the HFTs. We believe that the information ignored in this study could contain variables that would explain their order generation processes. Second, although this study only focused on the generation of limit orders by HFTs, it is also important to clarify the cancellation process for limit orders by HFTs and the generation of market orders. Third, we did not pay attention to profit and loss; however, practically, a key factor in an HFT strategy is the ability to make stable profits.
As the period of our data is very short, we did not observe any abnormal behavior in the market; however, we cannot deny the possibility that HFTs may overreact and result in serious synchronization during other periods or in other markets. Further studies of the relations among Hawkes parameters and the case of crashes are needed to prevent the excessive synchronization of biased orders of buy or sell. Our results are important since the model we derived in this paper provides a foundation for performing such studies through simulations. HFTs play a central role in providing liquidity to the market, and further detailed analyses of HFT strategies will contribute to the development of modern financial markets in general.

Author Contributions

Conceptualization, M.T.; formal analysis, H.W.; data curation, H.W.; investigation, H.W., H.T. and M.T.; methodology, H.W., H.T.; project administration, M.T.; resources, M.T.; software, H.W.; supervision, M.T.; writing—original draft, H.W.; writing—review and editing, H.T., M.T. All authors have read and agreed to submit this version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Restrictions apply to the availability of these data. Data were obtained from the EBS. We obtained permission for publication.

Acknowledgments

We appreciate EBS, NEX Group plc. for their provision of the EBS data.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Odean, T. Are investors reluctant to realize their losses? J. Financ. 1998, 53, 1775–1798. [Google Scholar] [CrossRef] [Green Version]
  2. Grinblatt, M.; Keloharju, M. The investment behavior and performance of various investor types: A study of Finland’s unique data set. J. Financ. Econ. 2000, 55, 43–67. [Google Scholar] [CrossRef]
  3. Sueshige, T.; Sornette, D.; Takayasu, H.; Takayasu, M. Classification of position management strategies at the order-book level and their influences on future market-price formation. PLoS ONE 2019, 14, e0220645. [Google Scholar] [CrossRef]
  4. Sueshige, T.; Kanazawa, K.; Takayasu, H.; Takayasu, M. Ecology of trading strategies in a forex market for limit and market orders. PLoS ONE 2018, 13, e0208332. [Google Scholar] [CrossRef] [PubMed]
  5. Kanazawa, K.; Sueshige, T.; Takayasu, H.; Takayasu, M. Kinetic theory for financial Brownian motion from microscopic dynamics. Phys. Rev. E 2018, 98, 052317. [Google Scholar] [CrossRef] [Green Version]
  6. Kanazawa, K.; Sueshige, T.; Takayasu, H.; Takayasu, M. Derivation of the Boltzmann equation for financial Brownian motion: Direct observation of the collective motion of high-frequency traders. Phys. Rev. Lett. 2018, 120, 138301. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  7. Biais, B.; Foucault, T.; Moinas, S. Equilibrium fast trading. J. Financ. Econ. 2015, 116, 292–313. [Google Scholar] [CrossRef] [Green Version]
  8. Carrion, A. Very fast money: High-frequency trading on the NASDAQ. J. Financ. Mark. 2013, 16, 680–711. [Google Scholar] [CrossRef]
  9. Schmidt, A.B. Ecology of the Modern Institutional Spot FX: The EBS Market in 2011. 2012. Available online: https://papers.ssrn.com/sol3/papers.cfm?abstract_id=1984070 (accessed on 7 January 2022).
  10. Gerig, A. High-Frequency Trading Synchronizes Prices in Financial Markets. 2015. Available online: https://papers.ssrn.com/sol3/papers.cfm?abstract_id=2173247 (accessed on 7 January 2022).
  11. Mukerji, P.; Chung, C.; Walsh, T. The impact of algorithmic trading in a simulated asset market. J. Risk Financ. Manag. 2019, 12, 68. [Google Scholar] [CrossRef] [Green Version]
  12. O’Hara, M. High frequency market microstructure. J. Financ. Econ. 2015, 116, 257–270. [Google Scholar] [CrossRef]
  13. Goldstein, M.A.; Kumar, P.; Graves, F.C. Computerized and high-frequency trading. Financ. Rev. 2014, 49, 177–202. [Google Scholar] [CrossRef] [Green Version]
  14. Jones, C.M. What Do We Know about High-Frequency Trading? Columbia Business School Research Paper 13-11. 2013. Available online: https://papers.ssrn.com/sol3/papers.cfm?abstract_id=2236201 (accessed on 7 January 2022).
  15. Menkveld, A.J. High frequency trading and the new market makers. J. Financ. Mark. 2013, 16, 712–740. [Google Scholar] [CrossRef] [Green Version]
  16. Easley, D.; De Prado, M.M.L.; O’Hara, M. The microstructure of the “flash crash”: Flow toxicity, liquidity crashes, and the probability of informed trading. J. Portf. Manag. 2011, 37, 118–128. [Google Scholar] [CrossRef] [Green Version]
  17. Easley, D.; de Prado, M.M.L.; O’Hara, M. VPIN and the flash crash: A rejoinder. J. Financ. Mark. 2014, 17, 47–52. [Google Scholar] [CrossRef]
  18. Andersen, T.G.; Bondarenko, O. VPIN and the flash crash. J. Financ. Mark. 2014, 17, 1–46. [Google Scholar] [CrossRef]
  19. Andersen, T.G.; Bondarenko, O. Reflecting on the VPIN dispute. J. Financ. Mark. 2014, 17, 53–64. [Google Scholar] [CrossRef]
  20. Andersen, T.G.; Bondarenko, O. Assessing measures of order flow toxicity and early warning signals for market turbulence. Rev. Financ. 2015, 19, 1–54. [Google Scholar] [CrossRef]
  21. D’Souza, C. Where Does Price Discovery Occur in FX Markets? 2007. Available online: https://papers.ssrn.com/sol3/papers.cfm?abstract_id=966446 (accessed on 7 January 2022).
  22. Gençay, R.; Gradojevic, N.; Olsen, R.; Selçuk, F. Informed traders’ arrival in foreign exchange markets: Does geography matter? Empir. Econ. 2015, 49, 1431–1462. [Google Scholar] [CrossRef] [Green Version]
  23. Gençay, R.; Gradojevic, N. Private information and its origins in an electronic foreign exchange market. Econ. Model. 2013, 33, 86–93. [Google Scholar] [CrossRef]
  24. Gradojevic, N.; Erdemlioglu, D.; Gençay, R. Informativeness of trade size in foreign exchange markets. Econ. Lett. 2017, 150, 27–33. [Google Scholar] [CrossRef]
  25. Elaut, G.; Frömmel, M.; Lampaert, K. Intraday momentum in FX markets: Disentangling informed trading from liquidity provision. J. Financ. Mark. 2018, 37, 35–51. [Google Scholar] [CrossRef]
  26. Hawkes, A.G. Spectra of some self-exciting and mutually exciting point processes. Biometrika 1971, 58, 83–90. [Google Scholar] [CrossRef]
  27. Engle, R.F.; Russell, J.R. Forecasting the frequency of changes in quoted foreign exchange prices with the autoregressive conditional duration model. J. Empir. Financ. 1997, 4, 187–212. [Google Scholar] [CrossRef]
  28. Takayasu, M.; Takayasu, H. Self-modulation processes and resulting generic 1/f fluctuations. Phys. A Stat. Mech. Its Appl. 2003, 324, 101–107. [Google Scholar] [CrossRef]
  29. Hawkes, A.G. Hawkes processes and their applications to finance: A review. Quant. Financ. 2018, 18, 193–198. [Google Scholar] [CrossRef]
  30. Bowsher, C. Modelling security market events in continuous time: Intensity based, multivariate point process models. J. Econom. 2007, 141, 876–912. [Google Scholar] [CrossRef] [Green Version]
  31. Filimonov, V.; Sornette, D. Quantifying reflexivity in financial markets: Toward a prediction of flash crashes. Phys. Rev. E 2012, 85, 056108. [Google Scholar] [CrossRef] [Green Version]
  32. Hardiman, S.J.; Bercot, N.; Bouchaud, J.P. Critical reflexivity in financial markets: A Hawkes process analysis. Eur. Phys. J. B 2013, 86, 442. [Google Scholar] [CrossRef] [Green Version]
  33. Hardiman, S.J.; Bouchaud, J.P. Branching-ratio approximation for the self-exciting Hawkes process. Phys. Rev. E 2014, 90, 062807. [Google Scholar] [CrossRef] [Green Version]
  34. Bacry, E.; Muzy, J.F. First- and Second-Order Statistics Characterization of Hawkes Processes and Non-Parametric Estimation. IEEE Trans. Inf. Theory 2016, 62, 2184–2202. [Google Scholar] [CrossRef]
  35. Bacry, E.; Jaisson, T.; Muzy, J.F. Estimation of slowly decreasing hawkes kernels: Application to high-frequency order book dynamics. Quant. Financ. 2016, 16, 1179–1201. [Google Scholar] [CrossRef]
  36. Achab, M.; Bacry, E.; Muzy, J.F.; Rambaldi, M. Analysis of order book flows using a non-parametric estimation of the branching ratio matrix. Quant. Financ. 2018, 18, 199–212. [Google Scholar] [CrossRef] [Green Version]
  37. Rambaldi, M.; Pennesi, P.; Lillo, F. Modeling foreign exchange market activity around macroeconomic news: Hawkes-process approach. Phys. Rev. E 2015, 91, 012819. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  38. Rambaldi, M.; Bacry, E.; Lillo, F. The role of volume in order book dynamics: A multivariate Hawkes process analysis. Quant. Financ. 2017, 17, 999–1020. [Google Scholar] [CrossRef]
  39. Rambaldi, M.; Filimonov, V.; Lillo, F. Detection of intensity bursts using Hawkes processes: An application to high-frequency financial data. Phys. Rev. E 2018, 97, 032318. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  40. Lu, X.; Abergel, F. High-dimensional Hawkes processes for limit order books: Modelling, empirical analysis and numerical calibration. Quant. Financ. 2018, 18, 249–264. [Google Scholar] [CrossRef] [Green Version]
  41. Fosset, A.; Bouchaud, J.P.; Benzaquen, M. Non-parametric estimation of quadratic Hawkes processes for order book events. Eur. J. Financ. 2021. [Google Scholar] [CrossRef]
  42. Rizoiu, M.A.; Lee, Y.; Mishra, S.; Xie, L. A tutorial on hawkes processes for events in social media. arXiv 2017, arXiv:1708.06401. [Google Scholar]
  43. Hawkes, A.G. Point Spectra of Some Mutually Exciting Point Processes. J. R. Stat. Soc. Ser. B Methodol. 1971, 33, 438–443. [Google Scholar] [CrossRef]
  44. Bacry, E.; Mastromatteo, I.; Muzy, J.F. Hawkes processes in finance. Mark. Microstruct. Liq. 2015, 1, 1550005. [Google Scholar] [CrossRef]
  45. Helmstetter, A.; Sornette, D. Importance of direct and indirect triggered seismicity in the ETAS model of seismicity. Geophys. Res. Lett. 2003, 30. [Google Scholar] [CrossRef] [Green Version]
  46. Helmstetter, A.; Sornette, D. Subcritical and supercritical regimes in epidemic models of earthquake aftershocks. J. Geophys. Res. Solid Earth 2002, 107, ESE-10. [Google Scholar] [CrossRef]
  47. Toke, I.M. An Introduction to Hawkes Processes with Applications to Finance. Lectures Notes from Ecole Centrale Paris, BNP Paribas Chair of Quantitative Finance. 2011, Volume 193. Available online: http://www.smallake.kr/wp-content/uploads/2015/01/HawkesCourseSlides.pdf (accessed on 7 January 2022).
  48. Kingma, D.P.; Ba, J. Adam: A method for stochastic optimization. arXiv 2014, arXiv:1412.6980. [Google Scholar]
  49. Ogata, Y. On Lewis’ simulation method for point processes. IEEE Trans. Inf. Theory 1981, 27, 23–31. [Google Scholar] [CrossRef] [Green Version]
  50. Ogata, Y. Statistical models for earthquake occurrences and residual analysis for point processes. J. Am. Stat. Assoc. 1988, 83, 9–27. [Google Scholar] [CrossRef]
  51. Lewis, P.W.; Shedler, G.S. Simulation of nonhomogeneous Poisson processes by thinning. Nav. Res. Logist. Q. 1979, 26, 403–413. [Google Scholar] [CrossRef]
  52. Ward, J.H., Jr. Hierarchical grouping to optimize an objective function. J. Am. Stat. Assoc. 1963, 58, 236–244. [Google Scholar] [CrossRef]
Figure 1. Schematic of trading in order book. In the EBS market, even a sell (buy) limit order becomes a hit buy (sell) if a buy (sell) limit order at the same price is already in the order book.
Figure 1. Schematic of trading in order book. In the EBS market, even a sell (buy) limit order becomes a hit buy (sell) if a buy (sell) limit order at the same price is already in the order book.
Entropy 24 00214 g001
Figure 2. (a) Average trading price per 10 min window; (b) daily number of orders for the 6 types of orders in unit of 10 5 .
Figure 2. (a) Average trading price per 10 min window; (b) daily number of orders for the 6 types of orders in unit of 10 5 .
Entropy 24 00214 g002
Figure 3. (a) Histogram of the minimum time intervals between orders for 134 HFTs individually; (b) hourly changes in the number of HFTs and non-HFTs participating in trading; and (c) hourly change of the percentage of limit orders provided by HFTs.
Figure 3. (a) Histogram of the minimum time intervals between orders for 134 HFTs individually; (b) hourly changes in the number of HFTs and non-HFTs participating in trading; and (c) hourly change of the percentage of limit orders provided by HFTs.
Entropy 24 00214 g003
Figure 4. Numbers of (a) sell limit orders and (b) buy limit orders per 10 min window for three HFTs (green: HFT with 4th highest order frequency; purple: HFT with 7th highest order frequency; yellow: HFT with 10th highest order frequency). (c) Numbers for six types of orders per 10 min window in the order book (red: sell limit order; blue: buy limit order; red dotted line: sell cancel; blue dotted line: buy cancel; orange: hit sell; sky blue: hit buy). The vertical axis of each figure shows the number of each type of order per 10 min window, and the horizontal axis shows the time from 0:00 to 18:00 on 6 June 2016.
Figure 4. Numbers of (a) sell limit orders and (b) buy limit orders per 10 min window for three HFTs (green: HFT with 4th highest order frequency; purple: HFT with 7th highest order frequency; yellow: HFT with 10th highest order frequency). (c) Numbers for six types of orders per 10 min window in the order book (red: sell limit order; blue: buy limit order; red dotted line: sell cancel; blue dotted line: buy cancel; orange: hit sell; sky blue: hit buy). The vertical axis of each figure shows the number of each type of order per 10 min window, and the horizontal axis shows the time from 0:00 to 18:00 on 6 June 2016.
Entropy 24 00214 g004
Figure 5. Schematic showing an example of the time evolution of the intensity function for a multivariate Hawkes process with the exponential kernel ϕ n , m ( t ) = α n , m e β n , m t .
Figure 5. Schematic showing an example of the time evolution of the intensity function for a multivariate Hawkes process with the exponential kernel ϕ n , m ( t ) = α n , m e β n , m t .
Entropy 24 00214 g005
Figure 6. Cumulative distributions of Δ t T S , i (a) and Δ τ ^ T S , i (b) for three typical HFTs (blue: HFT with 2nd highest order frequency; orange: HFT with 4th highest order frequency; and red: HFT with 7th highest order frequency).
Figure 6. Cumulative distributions of Δ t T S , i (a) and Δ τ ^ T S , i (b) for three typical HFTs (blue: HFT with 2nd highest order frequency; orange: HFT with 4th highest order frequency; and red: HFT with 7th highest order frequency).
Entropy 24 00214 g006
Figure 7. A histogram of D T S + T B K L values for all 134 HFTs. Out of 134 HFTs, 104 fell within the acceptable error threshold of 0.05, and the remaining 30 HFTs exceeded the threshold.
Figure 7. A histogram of D T S + T B K L values for all 134 HFTs. Out of 134 HFTs, 104 fell within the acceptable error threshold of 0.05, and the remaining 30 HFTs exceeded the threshold.
Entropy 24 00214 g007
Figure 8. Comparison between simulations and real data of { t T S , i } and { t T B , i } for HFT with D T S + T B K L = 0.0238 . The horizontal axis represents the time over a 24 h period, and the vertical axis represents the number of order occurrences per 10 min window.
Figure 8. Comparison between simulations and real data of { t T S , i } and { t T B , i } for HFT with D T S + T B K L = 0.0238 . The horizontal axis represents the time over a 24 h period, and the vertical axis represents the number of order occurrences per 10 min window.
Entropy 24 00214 g008
Figure 9. Comparison between simulations and real data of { t T S , i } and { t T B , i } for HFT with D T S + T B K L = 0.211 . The horizontal axis represents the time during a 24 h period, and the vertical axis represents the number of order occurrences per 10 min window.
Figure 9. Comparison between simulations and real data of { t T S , i } and { t T B , i } for HFT with D T S + T B K L = 0.211 . The horizontal axis represents the time during a 24 h period, and the vertical axis represents the number of order occurrences per 10 min window.
Entropy 24 00214 g009
Figure 10. Dendrogram based on the Ward method of clustering for 104 HFTs that were successfully modeled by the Hawkes process. The vertical axis represents the distance between the clusters, as defined in Equation (12), and the horizontal axis shows the labels of the HFTs according to the order frequency (red: Group A with 77 HFTs; blue: Group B with 12 HFTs; yellow: Group C with 15 HFTs).
Figure 10. Dendrogram based on the Ward method of clustering for 104 HFTs that were successfully modeled by the Hawkes process. The vertical axis represents the distance between the clusters, as defined in Equation (12), and the horizontal axis shows the labels of the HFTs according to the order frequency (red: Group A with 77 HFTs; blue: Group B with 12 HFTs; yellow: Group C with 15 HFTs).
Entropy 24 00214 g010
Figure 11. Percentile plot of normalized branching ratios (a) ρ ¯ T S , m and (b) ρ ¯ T B , m for 77 HFTs in group A. The vertical axis represents the normalized branching ratios by event m, and the horizontal axis represents element m M in both figures (top bar: 75th percentile; X symbol: median; bottom bar: 25th percentile; ◯ symbol: the mean).
Figure 11. Percentile plot of normalized branching ratios (a) ρ ¯ T S , m and (b) ρ ¯ T B , m for 77 HFTs in group A. The vertical axis represents the normalized branching ratios by event m, and the horizontal axis represents element m M in both figures (top bar: 75th percentile; X symbol: median; bottom bar: 25th percentile; ◯ symbol: the mean).
Entropy 24 00214 g011
Figure 12. Network graph of interaction between buy–sell limit orders of HFTs in Group A and all types of orders in the order book.
Figure 12. Network graph of interaction between buy–sell limit orders of HFTs in Group A and all types of orders in the order book.
Entropy 24 00214 g012
Figure 13. Percentile plot of normalized branching ratios (a) ρ ¯ T S , m and (b) ρ ¯ T B , m for 12 HFTs in Group B. The vertical and horizontal axes are the same as those in Figure 10 (top bar: 75th percentile; X symbol: median; bottom bar: 25th percentile; ◯ symbol: the mean).
Figure 13. Percentile plot of normalized branching ratios (a) ρ ¯ T S , m and (b) ρ ¯ T B , m for 12 HFTs in Group B. The vertical and horizontal axes are the same as those in Figure 10 (top bar: 75th percentile; X symbol: median; bottom bar: 25th percentile; ◯ symbol: the mean).
Entropy 24 00214 g013
Figure 14. Network graph of interaction between the buy–sell limit orders of HFTs in Group B and all types of orders in the order book.
Figure 14. Network graph of interaction between the buy–sell limit orders of HFTs in Group B and all types of orders in the order book.
Entropy 24 00214 g014
Figure 15. Percentile plot of normalized branching ratios (a) ρ ¯ T S , m and (b) ρ ¯ T B , m for 15 HFTs in Group C. The vertical and horizontal axes are the same as those in Figure 10 and Figure 11 (top bar: 75th percentile; X symbol: median; bottom bar: 25th percentile; ◯ symbol: the mean).
Figure 15. Percentile plot of normalized branching ratios (a) ρ ¯ T S , m and (b) ρ ¯ T B , m for 15 HFTs in Group C. The vertical and horizontal axes are the same as those in Figure 10 and Figure 11 (top bar: 75th percentile; X symbol: median; bottom bar: 25th percentile; ◯ symbol: the mean).
Entropy 24 00214 g015
Figure 16. Network graph of interaction between buy–sell limit orders of HFTs in Group C and all types of orders in the order book.
Figure 16. Network graph of interaction between buy–sell limit orders of HFTs in Group C and all types of orders in the order book.
Entropy 24 00214 g016
Table 1. Examples of raw data. Each order datum is tied to an anonymized trader ID.
Table 1. Examples of raw data. Each order datum is tied to an anonymized trader ID.
DateOrder TimeTrader IDOrder TypeUSD/JPYVolumeDeal Time
5 June 201621:00:12.946578Sell limit106.5151
5 June 201621:01:13.647HT6Buy cancel105.3902
5 June 201621:02:20.148JR1Buy limit105.405121:02:20.499
5 June 201621:02:20.499HSHSell market105.405121:02:20.499
5 June 201621:03:00.9507KPBid market106.4701
10 June 201620:59:20.148HT6Buy Limit107.405320:59:29.072
Table 2. Abbreviations for eight types of order events. Note the six types of orders in the order book do not include the orders of target HFTs. Therefore, { t i } differs for each HFT.
Table 2. Abbreviations for eight types of order events. Note the six types of orders in the order book do not include the orders of target HFTs. Therefore, { t i } differs for each HFT.
TS:Sell limit of the target HFT itselfTB:Buy limit of the target HFT itself
SL:Sell limit in order bookBL:Buy limit in order book
SC:Sell cancel in order bookBC:Buy cancel in order book
HS:Hit sellHB:Hit buy
Table 3. Mean values of the estimated time constant β ^ n , m 1 (s) for the HFTs in Group A.
Table 3. Mean values of the estimated time constant β ^ n , m 1 (s) for the HFTs in Group A.
mTSTBSLBLSCBCHSHB
n
TS0.1020.2950.0760.0870.0810.0690.1100.404
TB0.1180.0990.1140.0710.0660.0910.1480.111
Table 4. Sample means of 16 time constants, β n , m 1 (S), for HFTs in Group B.
Table 4. Sample means of 16 time constants, β n , m 1 (S), for HFTs in Group B.
mTSTBSLBLSCBCHSHB
n
TS0.0990.1000.6810.1090.1100.2850.1010.099
TB0.0990.1000.2590.3730.5580.1050.0990.101
Table 5. Sample means of 16 time constants, β n , m 1 (S), for HFTs in Group C.
Table 5. Sample means of 16 time constants, β n , m 1 (S), for HFTs in Group C.
mTSTBSLBLSCBCHSHB
n
TS0.13417.7810.2660.0990.0960.3590.1090.154
TB11.2010.0980.2690.2580.3550.0940.1320.103
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Watari, H.; Takayasu, H.; Takayasu, M. Analysis of Individual High-Frequency Traders’ Buy–Sell Order Strategy Based on Multivariate Hawkes Process. Entropy 2022, 24, 214. https://doi.org/10.3390/e24020214

AMA Style

Watari H, Takayasu H, Takayasu M. Analysis of Individual High-Frequency Traders’ Buy–Sell Order Strategy Based on Multivariate Hawkes Process. Entropy. 2022; 24(2):214. https://doi.org/10.3390/e24020214

Chicago/Turabian Style

Watari, Hiroki, Hideki Takayasu, and Misako Takayasu. 2022. "Analysis of Individual High-Frequency Traders’ Buy–Sell Order Strategy Based on Multivariate Hawkes Process" Entropy 24, no. 2: 214. https://doi.org/10.3390/e24020214

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop